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Abstract. We apply a semi-classical approach of handling waves as quasiparticle gas in a slowly varying flow - analogous to 
ray tracing - to calculate the Alfven wave transmission parameters, the resulting cross-helicity of the waves and the scattering- 
centre compression ratio, for cases where the shock thickness is large enough for the turbulent waves in the plasma to see the 
transition of the background flow parameters as smooth and slowly varying. For nonrelativistic shocks the wave transmission 
produces similar effects on the downstream turbulence and the scattering-centre compression ratio as does the transmission 
through a step shock: the downstream Alfven waves propagate predominantly towards the shock in the local plasma frame and, 
thus, the scattering-centre compression ratio is larger than the gas compression ratio. For thick relativistic shocks, however, we 
find qualitative differences with respect to the step-shock case: for low-Alfvenic-Mach-number shocks the downstream waves 
propagate predominantly away from the shock, and the scattering-centre compression ratio is lower than that of the gas. Thus, 
when taken into account, the Alfven wave transmission can decrease the efficiency of the first-order Fermi acceleration in a 
thick relativistic shock. 
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1. Introduction 

A common approximation in shock-related particle accelera- 
tion studies involving the first-order Fermi mechanism is to as- 
sume the scattering turbulence frozen-in to the plasma. This 
is very well justified in cases where the speed of the turbulent 
waves is negligibly small compared to the flow speed. For cases 
where the wave speed increases and becomes a notable fraction 
of the flow speed, however, the deviation of the scattering cen- 
tre speed from the speed of the background flow can have no- 
table effects on the particle acceleration scenario; it was noted 
already by Bell ( 1978 1 that the speed of the scattering centres 
is not necessarily that of the flow. Furthermore, the compres- 
sion ratio (of the scattering centres) is known to have signifi- 
cant effect on the power-law spectral index of particles acceler- 
ated in the shock for both steplike and thick shocks (for mod- 
ified shocks see Drury et al. 1982 and Drury 1983 for nonrel- 
ativistic analytical studies, and Virtanen & Vainio 2005a, here- 
after [VV05I, for numerical results for all speeds, and Keshet 
& Waxman 2005 for an analytical solution for step shocks). 
The turbulent waves are, of course, also affected by the shock, 
and thus one has to start from the transmission of turbulence 



through the shock to arrive at describing particle acceleration in 
a case where wave-propagation effects are taken into account. 



The effect of the transmission of Alfven waves through a 
step shock was solved at the nonrelativistic limit by Vainio & 
Schlickeiser ( 1998 1 and later generalised to include also rela- 
tivistic shocks by Vainio et al. (2003, hereafter referred to as 
IVVS03I and|5005l. They found that, regardless of the shock 
speed, waves initially in equipartition were propagating pre- 
dominantly antiparallel to the flow direction (i.e. backwards 
toward the shock in the downstream plasma frame) after being 
transmitted through a step shock with Alfvenic Mach number 
sufficiently low to yield non-negligible wave speeds compared 
to the speed of the plasma flow. They showed that this leads to 
an increased scattering-centre compression ratio, i.e., that the 
scattering-centre speed can undergo compression that is sig- 
nificantly higher than the compression of the flow itself. This 
effect thus causes the first-order Fermi acceleration to lead to 
accelerated particle energy spectral indices harder than ~ 2.2, 
the well-known outcome of parallel step shocks at the ultrarel- 
ativistic limit in case of isotropic particle scattering (e.g., Kirk 
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1 In Eq. (1) of VV05 there is a typo in the exponent in the denomi- 
nator of the last term: r 2 should be r 3 . 



& Duffy 1999, Keshet & Waxman 2005). This transmission 
analysis is, however, valid only for shocks of small thickness 
and waves with wavelengths much longer than the width of the 
shock transition. 
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In this paper we continue studying the transmission of 
Alfven waves through parallel shocks, but take an approach 
opposite to the step-shock approximation of VVS03] We cal- 
culate the transmission coefficients and their derivative param- 
eters (the average cross-helicity of the waves, and the compres- 
sion ratio of the scattering centres) for a continuous flow pro- 
file, approximating the transmission of waves through a shock 
with thickness sufficiently large for the waves to see the flow 
parameters varying slowly. The transmission coefficients for 
the nonrelativistic thick-shock case were calculated bv lVV05l 

In general, a single shock wave is not either a step-shock 
nor a thick shock: as the turbulence in the upstream medium 
has a broad-band spectrum, some waves will see the shock as 
a thick while others see it as a thin structure. In addition, some 
waves have wavelengths matching the shock thickness. Thus, 
the general theory of wave propagation can not rely on either 
the thick or thin shock approximation. However, the general 
theory is mathematically very cumbersome and we have here 
chosen to limit our discussion to waves with wavelengths much 
shorter than the shock thickness. Combining these results with 
our earlier ones ( V VS03 1 with the opposite assumption gives 
the first qualitative view on turbulence transmission through 
parallel shocks with finite thickness. 

2. Alfven wave transmission 

2.1. Theoretical background 

We adopt the idea of treating the waves as a quasiparticle gas. 
For monochromatic waves, the number density of the wave 
quanta, or the wave action density, is N' = E'/u>', where E' 
and u>' are the energy density and the angular frequency of the 
waves in the plasma frame. We will use the prime to denote the 
quantities measured in the plasma frame, while unprimed val- 
ues are measured in the shock frame. Below, we will consider 
a broad-band spectrum of waves and, thus, the evolution of the 
quasi-particle distribution function in the phase space (x, k). 

Our approach is equivalent to the use of a variational princi- 
ple in deriving the equation for the turbulent MHD wave-action 
density, pioneered by Dewar (1970 1 and extended to the rela- 
tivistic case by Achterberg ( 1983 1. 

In addition to the angular frequency u>, we will be needing 
the wavenumber k and frequency / = as well as their val- 
ues in the local plasma frame {to' , k! and /' = ^- respectively). 
Because the medium consists of a time-independent flow pro- 
file through the shock, the shock-frame frequencies are con- 
served during shock crossing and we can write the correspon- 
dence between the shock- and plasma-frame values as 



k = F(k' + V(J) = ±r(l ± VVaW/Va 



and 



/ = YW + Vk')/(2n) = ±T(V ± V A )f/V A , 



(1) 



(2) 



where the dispersion relation u>' = ±VpJ<.' is used. The Alfven 
speed Va is obtained from the proper Alfven speed (in units of 

c) as 



where p! = (p' + P')/n', n', p' and P' are the specific en- 
thalpy, the number density, the total energy density and the 
gas pressure, measured in the local plasma frame, and B' Q is 
the magnitude of the of the static large scale background mag- 
netic field. For easier comparison with the step-shock transmis- 
sion of VVS03 we will present results using the proper speed 
u = TV, where r = 1 / Vl - V 2 is the Lorentz factor, and the 
speeds are measured in units of c. Using subscripts 1 and 2 for 
far upstream and downstream values, respectively, we can in- 
troduce the Alfvenic Mach number of the shock M 2 — u\lu 2 hv 



We adopt the turbulence spectra in a power-law form 
/'*(&', x) = Iq ■ (k'Jk'f for k' > k' Q , 



(4) 



at wavelengths greater than k' Q , for waves flowing parallel (also 
referred to as forward waves, denoted by "+") and antiparallel 
(or backward, denoted by "-") to the flow direction. 

2.2. Equation for the wave flux 

Let us consider the propagation of an Alfven wave in a time- 
independent medium with spatial gradients only in the direc- 
tion of the background flow, aligned with the background mag- 
netic field and the x-axis. The Hamiltonian of the system is 



H(x,p) = h£l(x,p x /h), 



(5) 



where x and p = hk are the canonical coordinates and mo- 
menta, k is the wavevector, and 



a) — Q(x, k x ) — V g (x) k x 



(6) 



gives the dispersion relation of the Alfven wave. Here 
V g = (V + sVa)/(1 + sVVa) is the group speed of the wave and 
s — + 1 gives its propagation direction with respect to the back- 
ground plasma flow. Hamilton's equations of motion for the 
wave packet, thus, read 

. _ dH _ dQ. £ _ p x _ IdH _ dQ. 

dp x dk x ' H h dx dx' 

with all other equations of motion being trivial: y = z — p v — 
p z — because the Hamiltonian has no dependence on these 
phase space coordinates. Thus, the continuity equation for the 
phase space density, 

dJV 1 dN 1 

f(x > p) = ¥x^ = ¥¥x^k = ¥ N(x > k) > 

of the waves (or Liouville's equation, equivalently) yields 



(8) 



dt dx dk 
which, after the substitution of the equations of motion, reduces 
to 

— —(—m\-—(— 



N\ = 0. 



(9) 



We can integrate this equation over the perpendicular wave 
numbers k x and k z to obtain an equation for the number-density 
of waves over k x , 



m a = F a Va = B' J y/Anfi'n', 



(3) 



an rr 

(x,k x ) = = J J dk y dk z N(x,k), 
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as 



dn kx d tdQ. \ dim 



dt 



ox \ok x 



dk x \ dx 



0, 



(10) 



because Q does not depend on k y or k z . This equation holds 
regardless of the form of the distribution function with respect 
to the wavenumbers k y and k z . These are, of course, constants 
of motion for the waves in our system. 

From now on, we will assume that the wavevector is aligned 
with the x-axis (i.e., that k y — k z — 0) and simplify the nota- 
tion by writing k x = k. Thus, using the explicit form of the 
dispersion relation, we obtain in steady state 

d , v d (dv g 
dx^r ■dk[-dx- knk 



o. 



(ii) 



Consider, next, the wave-action density over frequency, 
n a = dN/(d 3 x d(L>) = n k /V g , and change to using the position 
x and the wave frequency a> - V g (x) k as the independent vari- 
ables instead of x and k. Thus, 



dx Va 



d_ 

dlo 



and we obtain 
Ml 



d_ 

(ho 

dVo 



~~dx~Vo 



d_ 



(Vg • V g n J 



= V, 



' d_ 

, dx^ 



In this equation a> appears only as a parameter numbering the 
modes, which can be all treated as monochromatic waves using 
the constancy of their flux 

V g n w = n k (x, to/VgixJ) (12) 
with respect to position at constant u. Thus, 

n k (x,k) = n k (-oo,ki), 

where k\(k,x) = V g (x)k/V gt i is the far-upstream wavenumber 
of the wave with wavenumber k at x. 

2.3. Wave transmission coefficients 

Next we calculate the transmission coefficients for Alfven 
waves with wavelengths much shorter than the shock thickness. 

Transmission coefficients are needed for solving the turbu- 
lence conditions at a given location in the shock. In contrast 
to the step shock case ( V VS03 1 where part of the waves are 
reflected at the shock and change their mode (from parallel 
to antiparallel, or vice versa), in the case of slowly changing 
medium all the waves are transmitted through the shock with- 
out reflection, and the turbulence at a given location can be 
obtained from equation 



I' ± (.k',x) = T' ± 2 (x)-I' 1 ± (k'), 



(13) 



where I'^ik') = I' ± (k',x 



-oo) are the turbulence spectra 



far upstream for parallel (+) and antiparallel (-) waves, and 
T'±{x) are the corresponding transmission coefficients, now to 
be solved. 

The co-moving intensity of the waves is related to the 
number-density of wave quanta by 



tui)'n' y ± 



4n 



V A 



V A k' tui> n\ 
~vfk 4tt 



ha>nf 



VjJT(l±VV A ) 4n r(V±V A ) 4n 

where we have made use of the invariance of under 
Lorentz transformations along the x axis. (This holds because 
of the Lorentz invariance of A^(jc, k) and of the perpendicular 
wavenumber element dk v dk z = dk[, dk' 7 .) Thus, we get 



nv±v A ) Il± 



l'*(x,k') = 



tkon* ri(Vi±V A i) 



4n 



Vai 



where T^Vi + Va.i)^ = cj = T(V + Va)^'- F° r power-law 
spectra, I'^ik') oc k'~ q , we finally obtain 



i> ± ( X ,k') = r l ± (k!)^- 

VA,\ 



ri(Vi ± Vai) 



T(V ± V A ) 



9+1 



(14) 



At the nonrelativistic limit in the far downstream, this re- 
duces to dVV05t 

19+1 



i^(k') = if(k') ■ t'I = if(k'y 



-9+1/2 



M+ 1 



M + r 1 ' 2 



(15) 



where r = V1/V2 is the gas compression ratio and M reduces 
to its nonrelativistic form, M = Vi/Va,i- 

3. Results and discussion 

From this point on, we deal only with plasma-frame quanti- 
ties, so we omit the primes in variable symbols and assume 
them being measured in the plasma rest frame, unless other- 
wise mentioned. The local plasma speeds V and the location 
are, of course, still in the shock frame. 

For a general view we first look at the case of vanishing up- 
stream cross-helicity. The ratio of forward-to-backward waves 
in the downstream is 

19+1 



F 



(Vi + Va,i)(V 2 - V A ,2) 



(Vi - V A ,i)(V 2 + V A ,2) 



(16) 



which, in contrast to transmission in step shocks, can be also 
greater than one (i.e., the downstream cross-helicity is posi- 
tive) if Va,2 < V^i/r. This, however, can occur only if the 
shock is relativistic (for nonrelativistic shock Va,2 = Va,i/ V^" 
and < 1 always). One immediately sees that, as a con- 

sequence, while the effect of the transmission in the nonrel- 
ativistic regime is very similar to that of both non- and fully 
relativistic step shocks, for relativistic shocks qualitative dif- 
ferences arise. 

In the following sections we study the cross-helicity of 
the transmitted waves in the downstream, and its development 
throughout the shock. We also calculate the scattering centre 
compression ratios in the downstream. 
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Fig. 1. The far downstream cross-helicity H C 2 as a function of 
Alfvenic Mach number M for different shock proper speeds 
Mi = V\Fi. M c = yfr is the critical Mach number, below which 
the downstream Alfven speed exceeds the local plasma speed, 
and the shock becomes non-evolutionary. The spectral index of 
the turbulence power-spectrum is q — 1.5 



3.1. Far downstream cross-helicity 

Here we calculate the normalised cross-helicity of the waves, 



H c (k,x) = 



I + (k, x) - r(k, x) 
I + (k, x) + I-(k, x)' 



(17) 



in the far downstream (// c> 2) from a given upstream state (H c i) 
as a function of different Alfvenic Mach numbers. 

For high-M shocks, where the effect of wave speeds is neg- 
ligible, the downstream cross-helicity approaches, of course, 
the cross helicity upstream as M — > oo. As the Alfvenic Mach 
number decreases and approaches the critical Mach number 
M c = yfr (i.e., as the downstream Alfven speed starts to ap- 
proach the downstream flow speed), the cross-helicity for non- 
relativistic shocks starts to decrease and approach - 1 (as shown 
bv lVV05t . as for step shocks of all speeds JVVS03t . 

For relativistic speeds and Mach number a few times the 
critical M c , however, nearly all of the waves are streaming par- 
allel to the flow and H c ^ — > +1 as u\ — > oo almost regardless of 
the upstream cross-helicity. For the lowest M > M c the cross- 
helicity drops to - 1 as the downstream Alfven speed becomes 
equal to the flow speed; this happens because the transmission 
coefficients of the backward waves, T\ goes to infinity. The 
physicality and removal of this mathematical singularity is dis- 
cussed in section l3~4l 

The lack of wave reflection causes the special cases of only 
one upstream wave field to keep also the downstream wave 



field at the same state (i.e., H { 



0,1 



±1 



:1) but the 



above-mentioned general behaviour is seen for all |// c ,il < 1. 
The downstream cross-helicity as a function of the Alfvenic 
Mach number is plotted in Fig^for non-to-highly relativistic 
shocks for vanishing upstream cross helicity. 



Fig. 2. Form of the flow profile (upper panel) and the local 
cross helicity across the shock when M/M c = 2 (lower panel). 
Shock proper speeds u\ range from 0.1 c to 10 c (lines from 
bottom to top, respectively). In the upper panel the flow speed 
is in arbitrary units, showing only the form of the transition 
from far upstream speed V\ to that of the downstream, Vj, and 
the spatial range across which the transition takes place for both 
the flow speed and the cross-helicities. See text for details. 



3.2. Cross-helicity as a function of location 

Next we want to study how the cross-helicity changes from the 
upstream value to that of the downstream. We use a hyperbolic 
tangent profile (Schneider & Kirk 1989 Virtanen & Vainio 
2005 1 to describe the flow speed as a function of location, 
and follow the development of the local mean cross-helicity. 
Since the actual shock thickness is of no importance for the 
spatial properties of the test-wave transmission, as long as it is 
large enough to allow for the "smoothly-varying-background" 
assumption, we use only one flow profile (and one thickness) to 
demonstrate the general behaviour. The form of the flow speed 
profile is shown in the top panel of Fig. |2] The unit of loca- 
tion x in that figure is arbitrary and is chosen so that the shock 
transition takes place within * 1 location unit. 

While for nonrelativistic shocks the cross-helicity changes 
monotonically from the upstream value to that calculated in 
Section Em for faster shocks the local cross-helicity H c (x) has 
non-monotonic behaviour: the cross-helicity forms a maximum 
near the "centre" of the flow profile, where the velocity gradient 
has its largest value, after which it decreases to the downstream 
value H C 2- 

An example of the non-monotonic cross-helicity is shown 
in Fig. [2] where the local cross-helicity H c {x) (here set equal to 
zero the far upstream) is plotted as a function of location x for 
several shocks with speeds ranging from non- to fully relativis- 
tic, but all having the Alfvenic Mach number M = 2M C . The 
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Fig. 3. Ratio of the scattering centre compression ratio r% to that 
of the gas r for different shock proper speeds u\ as functions of 
the Alvenic Mach number M. Lines are drawn for upstream 
cross-helicity H cl = -1 (group of the uppermost lines) and 
+ 1 (lower group). The spectral index of the turbulence power 
spectrum is q = 1.5 and M c = y7; see text for details. 

critical Mach number M c - yfr is calculated separately for 
each shock speed following the scheme described in VVS03 
The spectral index of the turbulence is q — 1 .5 for all cases. 

3.3. Decreased compressio ratio 

While the assumption of small-amplitude waves leaves the 
compression of the plasma itself unaffected, the compression 
felt by the scattering particles, on the other hand, can be 
changed remarkably. This is due to the fact that the particles are 
scattered by the turbulent waves, and if the waves have speed 
non-neglibigle compared to the speed of the flow (i.e. waves are 
not assumed to be frozen-in to the plasma), the speed profile of 
the scattering centres can be different from that of the plasma. 
Thus, also the effective compression ratio felt by the particles, 
can differ from the gas compression ratio r. Using the local 
cross-helicity H c we can express the average wave speed as 



V k (x) 



V(x) + H c (x)V A (x) 
1 + H c (x)V(x)V A (xy 



(18) 



Now the scattering-centre compression ratio is the calculated 
using the average wave speeds (V^i and V^i)'- ru = Vt,i/V*,2- 

As noted earlier, if the cross-helicity in the upstream equals 
+ 1 it remains unchanged throughout the shock. So if there are 
e.g. only forward waves in the upstream, the upstream scat- 
tering centre speed is simply the (relativistic) sum of the flow 
speed and the local Alfven speed; and likewise in the down- 
stream. However, because of the compression and drop in the 
flow speed at the shock, the ratio of the Alfven and the flow 
speed in the downstream is larger than in the upstream. Also 
its effect on the resulting shock-frame scattering-centre speed 
is larger when the underlying flow speed is not close to the 



Fig. 4. Same as Fig. [3 but for H cl = -0. 



speed of light. For only forward waves in the upstream this re- 
sults in ru < r, and for only backward waves in the upstream 
Tk > r. When the Alfven speed drops and becomes negligible 
compared to local flow speeds (i.e. M — > oo), — > r. This is 
illustrated in Fig. [3] 

For vanishing upstream cross-helicity, the behaviour in a 
nonrelativistic shock was examined by V V05 and it was shown 
to be very similar to that in step shocks. Now for thick relativis- 
tic shock, because of the positive downstream cross-helicity, 
the average shock-frame wave speed in the downstream is 
higher than the local flow speed, and, thus, the compression 
ratio of the scattering centres is lower than that of the gas. 

Resulting scattering-centre compression ratios (scaled to 
the gas compression ratio) are plotted for upstream cross- 
helicities H c j = +1 in Fig- ar, d f° r H c l = in Fig.0] The 
nonrelativistic (u\ = 0.1) case is the same as in VV05 and, as 
for step shocks, the compression ratio increases also for higher 
speeds. However, on the contrary to what lVV05l expected. for 
the fastest low-M shocks the scattering-centre compression ra- 
tio does not increase steadily as M — > M c , but falls below the 
gas compression ratio when the shock speed becomes relativis- 
tic. 

3.4. Limitations for the analysis 



In earlier studies by Vainio & Schlickeiser (1998) for non- 
relativistic shocks and IVVS0 3 for relativistic ones, the trans- 
mission properties were solved for waves of wavelengths es- 
sentially longer than the length scales of the shock. Here we 
have studied the opposite of "short" waves flowing through a 
shock thicker than the wave length scales. Although these two 
transmission schemes work differently and in different parts of 
parameter space, it would be physically meaningful to apply 
both schemes in the same shock: the analysis o flVVS0 3 for the 
part of the turbulence spectrum for which the wavelengths are 
longer than the width of the shock transition, and the method 
described in this paper for waves for which the plasma parame- 
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Fig. 5. Sketch of the scattering-centre compression ratio as a 
function of wavenumber (left) and the resulting accelerated 
particle spectrum (right) in a thick relativistic shock. L denotes 
the shock thickness. 

ters vary slowly throughout the transition. Thus, the scattering- 
centre compression ratio is a function of wavenumber, and the 
resulting accelerated particle spectrum is not a simple power- 
law in energy (Fig. O- 

As the study here assumed no back-reaction of the waves 
to the shock structure or dynamics, the approach is limited to 
small-amplitude waves. The test-wave approach, on the other 
hand, leads to mathematical singularities at M — > M c — -\fr; at 
this limit the antiparallel waves are amplified infinitely. In the 
analysis this limitation affects through the simplified calcula- 
tion of the gas compression ratio. In more detailed further anal- 
ysis the effect of the upstream wave parameters to the calcula- 
tion of the gas compression ratio should be taken into account 
in order to remove the (non-physical) singularities, as demon- 
strated by Vainio & Schlickeiser ( 1999 1 for nonrelativistic step 
shocks. 

Also the inclusion of accelerating particles will have effects 
on the waves; in this treatise the effect of the particles on the 
waves have been omitted. However, the resonant wave-particle 
interactions could have a damping effect on the turbulence, thus 
rapidly changing the transmitted wave distributions, especially 
when moving downstream away from the shock. Similarily, our 
analysis also neglects the wave-wave interactions in the down- 
stream medium (Vainio & Spanier 2005 1. 

Finally we merely note that in the present study we have 
confined the analysis to strictly parallel shocks. For oblique 
magnetic field alignments the treatment of the turbulence trans- 
mission becomes very complicated and is beyond the scope of 
this study. 



(i.e. that of the scattering centres) at a thick relativistic low- 
M shock was shown to be weaker than the compression of the 
background flow. For slower (from non- to mildly relativistic) 
shocks the scattering centre compression ratio was shown to 
increase, like in step shocks. 
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4. Conclusions 

We have calculated the transmission coefficients for Alvfen 
(test) waves with short wavelengths through a modified shock 
with background plasma properties changing slowly compared 
to the wave length scales. 

While the transmission trough a step shock was found sim- 
ilar throughout the whole shock speed range ( V VS03 1, and, ad- 
ditionally, very similar behaviour was observed also for thick 
shocks in the nonrelativistic limit ( V V05 1, qualitative and sig- 
nificant differences emerge for thick shocks when the speed 
increases. As a consequence, from the point of view of the par- 
ticle acceleration studies, the compression felt by the particles 



